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A measurement of the top quark mass using events with one charged lepton, missing transverse 
energy, and jets in the final state, collected by the DO detector from pp collisions at = 1.96 
TeV at the Fermilab Tevatron collider, is presented. A constrained fit is used to fully reconstruct 

the kinematics of the events. For every event a top quark mass likelihood is calculated taking 
into account all possible jet assignments and the probability that an event is signal or background. 
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Lifetime-based identification of b jets is employed to enhance the separation between tt signal and 
background from other physics processes and to improve the assignment of the observed jets to the 
quarks in the tt hypothesis. We extract a multiplicative jet energy scale factor JES in-situ, greatly 
reducing the systematic effect related to the jet energy measurement. In a data sample with an 
integrated luminosity of 425 pb~^, we observe 230 candidate events, with an estimated background 
of 123 events, and measure mt = 173.7±4.4 (stat + JES)]'^2 J (syst) GeV. This result represents the 
first application of the Ideogram technique to the measurement of the top quark mass in lepton+jets 
events. 

PACS numbers: 14.65.Ha, 12.15.Ff 



I. INTRODUCTION 

The top quark is by far the heaviest known fermion. 
Its discovery in 1995 [T] confirmed the structure of the 
standard model, while its strikingly large mass compared 
to other fermions highlights remaining open questions re- 
lated to the large range of quark and lepton masses and 
the precise mechanism of electroweak symmetry breaking 
that explains fermion masses in the theory. Within the 
framework of the standard model, the top quark mass is 
related to the Higgs boson mass and the W boson mass 
through radiative corrections. A precise measurement 
of the top quark mass helps to constrain the standard 
model and to predict the mass of the Higgs boson 0. 
At the same time it provides a challenge to the standard 
model with increased precision and distinguishes possible 
extensions of it. 

At the Fermilab Tevatron Collider, which collides pro- 
tons and antiprotons with a center-of-mass energy of 
1.96 TeV, the top quark is predominantly produced in tt 
pairs through qq annihilation 85%) and gluon-gluon 
fusion. In the framework of the standard model, the top 
quark decays almost exclusively to a 6 quark and a W 
boson. Thus the final state topology of a tt event is 
determined by the decay modes of the two W bosons. 
The analysis presented in this paper uses the lepton-|-jets 
(£-|-jets) channel, where one W boson decays hadronically 
and the other W boson decays leptonically to a muon or 
an electron and the corresponding (anti)neutrino. Tau 
leptons are not explicitly reconstructed in the analysis. 
Throughout this paper, charge conjugate modes are im- 
plicitly included. 

The £-|-jets topology combines a sizable branching frac- 
tion with a striking signature of the isolated energetic lep- 
ton and large missing transverse energy from the escaping 
neutrino. The background from M^-t-jets and QCD mul- 
tijet events is manageable. This means that the £+jets 
channel is particularly suited for studies of top quark 
properties, and it has provided the most precise mea- 
surements of the top quark mass to date 0,13, B B 0| ■ 

Nonetheless, serious challenges exist. A direct mea- 
surement of the top quark mass requires that the kine- 
matics of the event are fully reconstructed, including the 
momentum of the neutrino. The signal events need to be 
separated from backgrounds in a manner that does not 
bias the mass measurement. Furthermore, with four jets 
in the final state, the assignment of jets to the original 



top quark decay products gives a twelve-fold ambiguity. 
Finally, a proper calibration of the jet energy scale is cru- 
cial. In previous measurements, this was the dominant 
systematic uncertainty. 

Early measurements of the top quark mass 3] (and 
some recent analyses Q) used a constrained fit to re- 
construct the kinematics of the event, choosing one jet 
assignment based on the quality of the fit. A distribu- 
tion of some event variable strongly correlated to the top 
mass, typically the fitted top mass, was plotted for data 
events. This distribution was then compared to distri- 
butions based on Monte Carlo simulation generated for 
different top quark masses to determine the value of the 
top quark mass that best agrees with the data. In the 
case of the DO analysis Q, a multi-variate discriminant 
that separates signal from background was also used in 
a two-dimensional likelihood fit to the Monte Carlo ref- 
erence distributions. However, in these analyses, only a 
certain amount of information per event is used in the 
final fit. 

The DO Matrix Element analysis Q demonstrated for 
the first time that the statistical precision of the measure- 
ment can be greatly enhanced by constructing event-by- 
event likelihoods that reflect the full ambiguity of the 
events. A dramatic improvement was achieved, albeit at 
the cost of computationally intensive methods. 

The analysis presented in this paper uses the Ideogram 
technique. This method is based on a constrained kine- 
matic fit and strives to obtain a similar improvement in 
statistical precision as the Matrix Element analysis with 
minimal additional computation. The constrained fit is 
used to determine the kinematics of the events and to 
improve their reconstruction beyond the detector resolu- 
tion. A top quark mass likelihood is derived for every 
event including all possible assignments of jets to quarks 
in the tt hypothesis, and taking into account the possi- 
bility that the event is background. The top quark mass 
is extracted through a combined likelihood fit including 
all events. This approach is very similar to the Ideogram 
technique used by the DELPHI experiment to measure 
the W boson mass at the CERN LEP collider 8]. Also 
there the different possible jet permutations lead to an 
ambiguity in the mass fit which is reflected in the event 
likelihood as the sum of Gaussian resolution functions. 
The similarity with the ideogram plots used by the Par- 
ticle Data Group Q to visualize a set of measurements 
is what gave the method its name. This is the first time 
the method is used to determine the top quark mass in 
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the £+jets channel. Recently, it has also been applied to 
the all-hadronic decay channel [To| . 

The free parameters in the fit are the top quark mass, 
the tt signal fraction in the sample, and an overall jet en- 
ergy scale (JES) factor. Including the JES factor as a free 
parameter in the fit greatly reduces the systematic un- 
certainty related to the jet energy scale calibration [1, H]- 
We employ the tagging of 6 jets, i.e., 6 tagging, to enhance 
the separation between signal and background from other 
physics processes. The b tags also help to better distin- 
guish between correct and wrong jet assignments in the 
likelihood. Events with and without b tags are included 
in the overall likelihood fit. 

This paper is organized as follows: Sections HIl and Hill 
describe the DO Run II detector and the event reconstruc- 
tion, respectively. Sections HV] to |Vl] describe the data 
and simulation samples used and outline the event selec- 
tion. In Sec. IVIIl the sample composition is estimated 
using topological and b tagging information. Section|VlTT] 
describes in detail the calculation of the Ideogram like- 
lihood and the Monte Carlo calibration procedure. The 
method is applied to data in Sec. IIXI and the systematic 
uncertainties are discussed in Sec.|Xl Section IXll presents 
a cross-check of the JES calibration, followed by the con- 
clusion in Sec. IXIII 



followed by fine and coarse hadronic sections with mod- 
ules assembled in a projective geometry to the interaction 
region. In addition to the preshower detectors, scintilla- 
tors between the CC and EC cryostats provide sampling 
of developing showers for 1.1 < \ri\ < 1.4. 

A muon system (l3| resides beyond the calorimetry and 
consists of a layer of tracking detectors and scintillation 
trigger counters before 1.8 T iron toroids, followed by 
two similar layers after the toroids. Tracking for jr/j < 1 
relies on 10 cm wide drift tubes [l2|, while 1 cm mini-drift 
tubes are used for 1 < j?;! < 2. 

Trigger and data acquisition systems are designed to 
accommodate the high luminosities of Run II. Based on 
preliminary information from tracking, calorimetry, and 
muon systems, the output of the first level of the trigger 
is used to limit the rate for accepted events to approx- 
imately 2 kHz. At the next trigger stage, with more 
refined information, the rate is reduced further to about 
1 kHz. These first two levels of triggering rely mainly on 
hardware and firmware. The third and final level of the 
trigger, with access to all of the event information, uses 
software algorithms and a computing farm, and reduces 
the output rate to about 50 Hz, which is written to tape. 

III. EVENT RECONSTRUCTION 



II. THE DO DETECTOR 

The Run II of the Fermilab Tevatron collider started 
in 2001 after substantial detector upgrades following the 
first Tevatron coUider run in 1992-1996. The DO Run II 
detector [ll| consists of a magnetic central tracking sys- 
tem, composed of a silicon micro-strip tracker (SMT) and 
a central fiber tracker (CFT), both located within a 2 T 
superconducting solenoidal magnet. The SMT has ap- 
proximately 800, 000 individual strips, with typical pitch 
of 50 — 80 /im, and a design optimized for tracking and 
vertexing capabilities at pseudorapidities of \r]\ < 2.5. 
The system has a six-barrel longitudinal structure, each 
with a set of four layers arranged axially around the beam 
pipe, and interspersed with 16 radial disks. The CFT has 
eight thin coaxial barrels, each supporting two doublets 
of overlapping scintillating fibers of 0.835 mm diameter, 
one doublet being parallel to the beam axis, and the other 
alternating by ±3° relative to the axis. Light signals are 
transferred via clear fibers to solid-state photon counters 
(VLPCs) that have « 80% quantum efficiency. 

Central and forward preshower detectors located just 
outside of the superconducting coil (in front of the 
calorimetry) are constructed of several layers of extruded 
triangular scintillator strips that are read out using 
wavelength-shifting fibers and VLPCs. The next layer of 
detection involves three liquid- argon/uranium calorime- 
ters: a central section (CC) covering approximately |?7| < 
1.1, and two end calorimeters (EC) that extend coverage 
to \?]\ « 4.2, all housed in separate cryostats [l2|. The 
calorimeter consists of an electromagnetic (EM) section 



This section summarizes the offiine event reconstruc- 
tion. We use a right-handed Cartesian coordinate sys- 
tem with the z axis defined by the direction of the pro- 
ton beam, the y axis pointing vertically upwards and the 
X axis pointing out from the center of the accelerator 
ring. The origin is at the center of the detector. The 
polar angle is defined with respect to the positive z 
axis and (j) is the azimuthal angle from the x axis in the 
transverse xy plane. The pseudorapidity r] is defined as 
T] = -ln(tan(6l/2)). 

A. Tracks and event vertex 

Tracks are reconstructed from the hit information in 
the SMT and CFT. A Kalman filter [14] is used to fit 
track candidates found by a road-based algorithm or 
a technique searching for clusters of track parameters 
formed by tracker hits. Using a vertex search procedure 
[l5| . a list of reconstructed primary vertices is returned. 
The primary event vertex for the tt reconstruction is cho- 
sen from this list based on the pt spectrum of the tracks 
associated with a given vertex. Only vertices with at 
least three tracks associated with them are considered. 



B. Electrons 

We reconstruct electrons using information from the 
calorimeter and the central tracker. Clusters of EM 
calorimeter cells (EM clusters) are built with a simple 
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cone algorith m using seeds o f Et > 1-5 GeV and ra- 
dius AR = v/(A?7)2 + (A0)2 = 0.2. An "extra-loose" 
electron is defined as an EM cluster with 90% of its en- 
ergy from the EM part of the calorimeter and isolated 
from hadronic energy depositions. Its longitudinal and 
transverse energy profiles have to be consistent with ex- 
pectations from simulation. In addition, the electrons 
used in the final event selection are required to match 
a track reconstructed in the central tracker and to pass 
an electron likelihood cut. The likelihood is built from 
seven variables containing tracking and calorimeter in- 
formation and is optimized to discriminate between true 
electrons and background. 



C. Muons 

Muons are reconstructed from the information in the 
muon system and the central tracker. We require a muon 
candidate to have hits in the muon detectors both in- 
side and outside the toroid. The timing information of 
the scintillator hits has to be consistent with that of a 
particle produced in a pp collision, thus rejecting cosmic 
muons. The muon candidate track is then extrapolated 
to the point of closest approach to the beam line, and 
matched to a track from the central tracking system us- 
ing a global track fit. Muons must not be surrounded by 
activity in the tracker or calorimeter and are required to 
be separated from reconstructed jets by AR > 0.5. 



D. Jet reconstruction and energy scale 

Jets are reconstructed from the calorimeter informa- 
tion using a cone algorithm [l6[ with radius AR = 0.5. 
Only calorimeter cells with signal larger than 4(t above 
the average noise and adjacent cells with signal at least 
2(7 above the noise are used. The jets are required to be 
confirmed by independent calorimeter trigger informa- 
tion and must be separated from an extra-loose electron 
by AR > 0.5. The reconstructed jet energies iSjeco 
corrected for an energy offset Eog, energy response i?cai, 
and out-of-cone showering Ccone, according to: 



rcco 



^cal ' Co 



(1) 



The offset correction is determined from events taken 
with a zero bias trigger during physics data taking and 
accounts for noise, multiple interactions, and energy pile- 
up. The response correction is derived from a high statis- 
tics jet+photon sample by looking at the px imbalance 
in these events. The photon energy scale is assumed to 
be equivalent to the well-known electron energy scale as 
calibrated from Z ee events. The showering correc- 
tion accounts for energy that particles inside the jet cone 
deposit outside the cone during the hadronic showering 
process. Transverse jet energy profiles are studied to de- 
termine this correction. Jets containing a muon within 



the jet cone are further corrected for the momentum car- 
ried by the muon and the associated neutrino. Since the 
method to extract the top quark mass is calibrated with 
respect to the Monte Carlo simulation, it is important to 
determine the relative jet energy scale S between data 
and the Monte Carlo simulation, 



S = 
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S is parameterized as a function of photon pT for several 
bins in (ptiV) space and is found to be flat within its un- 
certainties. No corrections from this source are therefore 
applied. Effects of a potential pt dependence are taken 
into account as a systematic uncertainty. For the overall 
jet energy scale, a uniform factor, JES, is introduced as 
a free parameter in the analysis. This factor is fitted in 
situ, simultaneously with the top quark mass in data by 
using information from the invariant mass of the hadron- 
ically decaying W bosons. For every event, this mass is 
constrained in the kinematic fit to be equal to the known 
value of the W boson mass ^] . The of the kinematic 
fit reflects the compatibility of the reconstructed jet en- 
ergies with this constraint. The likelihood is sensitive 
to the JES parameter through the x^- The overall fit 
will give the maximum likelihood for the value of JES 
which optimizes, on average, the compatibility between 
the reconstructed and fitted jet energies. 

Apart from the W boson mass information, no con- 
straint on the overall energy scale is used in the top quark 
mass fit. The jet energy scale measured in situ is con- 
sistent with the result obtained from photon+jet studies 
(SeclU). 

The analysis is calibrated such that in pseudo- 
experiments with Monte Carlo events the average fitted 
JES value is equal to one. A fitted value JES < 1 means 
that the jet energies in the sample considered are under- 
estimated with respect to the reference Monte Carlo scale 
described above {JES < 1 is equivalent to 5 < when 
fitting the data sample). 



E. Missing transverse energy 

We identify neutrinos indirectly from the energy imbal- 
ance in the event. The imbalance is reconstructed from 
the vector sum of the transverse energies in the calorime- 
ter cells and the reconstructed muons. Energies from the 
cells in the coarse hadronic portion of the calorimeter are 
only added if associated with a reconstructed jet. The 
missing transverse energy, ^t , is corrected for the en- 
ergy scale calibration of jets and electrons. 



F. 6-jet identification 

We identify b jets using a lifetime tagging algorithm 
(secondary vertex tagger, SVT) based on the explicit re- 
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construction of a secondary vertex from the decay of a 
&- flavored hadron • We call dca the distance of closest 
approach between a track and the beam line, with a{dca) 
being the uncertainty on dca. After the reconstruction 
of the primary event vertex, we consider tracks with 
dca/ a[dca) > 3.5 for the reconstruction of additional 
(secondary) vertices. For a reconstructed secondary ver- 
tex, the transverse decay length L^y with respect to the 
primary event vertex is computed. A jet is tagged as a 6 
jet if a secondary vertex is reconstructed within AR < 0.5 
of the jet with L^y/ a^L^y) > 7.0, where a{Lxy) is the un- 
certainty on Lxy The 6-jet tagging rate et is measured 
in data using information from an independent ^-tagging 
analysis that looks for the presence of a muon in the jet 
cone. 

Light quark jets can also be tagged when a fake 
secondary vertex is reconstructed due to track mis- 
measurements and random overlaps of tracks. This light 
jet tagging rate e; is estimated from the rate of secondary 
vertices with Lxy/(j{Lxy) < —7.0 in a data sample with 
predominantly light quark jets. Negative values of L^y 
occur if the secondary vertex is on the opposite side of 
the event vertex with respect to the jet and are a sign 
of mis-measurement and resolution effects. Misrecon- 
structed vertices with negative and positive values of L^y 
are expected to occur at the same rate. Corrections for 
the contamination with heavy flavor and the presence of 
long lived particles are applied as determined from Monte 
Carlo simulation. The 6-jet and light-jet tagging rates are 
measured in data and are parametrized as a function of 
jet transverse momentum and pseudorapidity [l5l |. 



muon-f jets channel. Similarly, in the electron-|-jets chan- 
nel the candidate electron must not be matched to a track 
or fail to pass the electron likelihood cut (Sec. IIIIBjl . 

V. SIMULATION 

Monte Carlo event generators are used to create large 
samples of simulated signal and background events. 
These samples are used for the calibration of the cen- 
tral mass value and the estimate of the uncertainty. We 
use ALPGEN 1.3 18] to generate signal and M^-|-jets back- 
ground events. The underlying event and hadronization 
is simulated using pythia 6.2 [l^. Signal tt events are 
generated at nine mass points with masses ranging from 
150 GeV to 200 GeV. The factorization and renormal- 
ization scales are set to Q = mt for the tt simulation 
and Q2 = _|_ ^(pj^')2 for W+\ets. All events are 
passed through a full GEANT-based [20| DO detector sim- 
ulation and reconstructed with the same software as the 
collider data. Events are accepted according to the prob- 
ability that a simulated event would pass the trigger re- 
quirements. This probability is typically between 0.9 and 
1.0. The same object and event selections as for the data 
samples are applied. The simulation chain is tuned to re- 
produce resolutions of reconstructed objects seen in the 
collider data. 



VI. KINEMATIC FIT AND FINAL EVENT 
SELECTION 



IV. DATA SAMPLES 

This paper describes the analysis of data collected be- 
tween April 2002 and August 2004, corresponding to an 
integrated luminosity of approximately 425 pb~^. For 
this analysis, the data sample was selected by triggering 
on a lepton and at least one additional jet in the events. 
The specific t rigg er requirements are described in more 
detail in Ref. p^ . 

The event selection requires an isolated lepton of trans- 
verse momentum pT > 20 GeV, with a pseudorapidity 
\ri\ < 1.1 for electrons and \ri\ < 2 for muons. Missing 
transverse energy fJx > 20 GeV is required as well as 
four or more jets with > 20 GeV and jryj < 2.5. A 
A(j) cut between fJx and lepton momentum is imposed 
to exclude events where the transverse energy imbalance 
is caused by a poor measurement of the lepton energy. 
The position of the event vertex along the beam direction 
has to be within 60 cm of the center of the detector. We 
select 246 candidate events. 

A QCD multijet background sample is also extracted 
from data by reversing the final lepton quality require- 
ment. Leaving all other event selection cuts unchanged, 
the candidate isolated muon must fail to be isolated from 
activity in the tracker or calorimeter (Sec. IIII C|l in the 



The kinematics of the events, including the undetected 
neutrino from the W boson decay, are reconstructed us- 
ing the same kinematic constrained fit that was devel- 
oped for the Run I analysis ji] . The resolutions of muons, 
electrons and jets were updated for Run II 0, IMlHSl- 

In events with more than four jets, only the four jets 
with highest px arc considered as possible candidates to 
be a light quark or b quark in the ti hypothesis used in 
the constrained fit. 

All twelve possible assignments of jets to quarks are 
considered. As a starting point for the kinematic fit, the 
unmeasured component of the neutrino momentum par- 
allel to the beam, p^, is chosen such that the two top 
quarks are assigned equal mass. This yields a quadratic 
equation for p^. We use both solutions as input to the fit 
yielding twenty-four fit results per event. Depending on 
the event kinematics and resolution effects, the discrimi- 
nant of the quadratic equation may be negative, in which 
case the discriminant is forced to be zero. Thus one or 
two solutions are always obtained. If only one solution 
is available, we include the same fit result twice in the 
likelihood. 

For the kinematic fit, we relate the reconstructed jet 
energy to the unfragmented parton energy. To this end, a 
jet-parton energy mapping is applied, which is the same 
in data and MC simulation. The corrections depend on 
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the flavor (b quark or light quark) of the parent quark 
and therefore depend on the jet-to-parton assignment 
used. To derive the mapping functions, we use MC events 
where the jets are unambiguously matched to the partons 
of the tt decay and compare the jet energy to the MC gen- 
erated parton energy. The jet-parton mapping functions 
contain the JES parameter as a uniform multiplicative 
factor. 

The kinematic fit is performed by minimizing a 
subject to the kinematic constraints: m{t — > Ivb) — 
mit qqb), m{£iy] = M^y, and 'm{qq) = M^. We use 
M\Y = 80.4 GeV 9]. The minimization algorithm uses 
the method of Lagrange multipliers; the nonlinear con- 
straint equations are solved using an iterative technique. 
From the fit for each jet/neutrino solution i, we extract 
the mass rrii , the estimated uncertainty on the fitted mass 
CTj, and the goodness of fit Xi- The fit is repeated for dif- 
ferent values of JES. The JES parameter is varied in 
steps of 3% in an interval of ±15% around unity. Only 
jet combinations for which the fit converges at all values 
of JES are used. This requirement is needed to prevent 
discontinuities as a function of JES in the event likeli- 
hood. The fitted mass mi(JES), estimated uncertainty 
ai{JES), and goodness of fit xfiJES) all depend on the 
JES parameter. In the following this dependence is not 
shown explicitly, to improve readability. 

The final selection requirement is that at least one 
jet/neutrino solution yields < 10 for the kinematic 
fit with JES = 1. This cut reduces the number of events 
from 120 to 116 in the electron+jets channel and from 
126 to 114 in the muon-|-jets channel. Most of the events 
removed by this cut are background events or badly re- 
constructed tt events that do not satisfy the tt fit hypoth- 
esis and do not carry useful information about the top 
quark mass. The algorithmic efficiency of the kinematic 
fit is excellent, as listed in Table [J 



TABLE I: The numbers of events and efficiencies for the elec- 
tron-|-jets (e) and muon-|-jets (/i) channel having at least one 
jet combination for which the fit converges at JES = 1, with- 
out and with the requirement on the maximum value of the 
X^- Each column shows the relative efficiency with respect to 
the previous column. 

Convergence of the kinematic fit: 





before 


converges, 


x' < 10 


converges 




fit 


JES = 1 


JES = 1 


aU JES* 


tt e 


9452 


100.0% 


97.7% 


100.0% 


M 


9265 


99.8% 


94.4% 


100.0% 


VK-l-jets e 


5163 


100.0% 


94.2% 


100.0% 


M 


5820 


99.7% 


89.9% 


100.0% 


data e 


120 


100.0% 


97.0% 


100.0% 


M 


126 


100.0% 


91.0% 


100.0% 



• for aU values of JES in the fit range 0.85 < JES < 1.15 



VII. SAMPLE COMPOSITION 

In order to obtain a good separation between tt sig- 
nal and background events (mainly VF-l-jets), a likelihood 
discriminant based on the 'low-bias' topological discrim- 
inant Dlb, developed in Run I p], is used. The -Dlb dis- 
criminant was designed to have minimal correlation with 
the top mass and is based on the following four topolog- 
ical variables: , aplanarity, H^2i ^'Hd. K^^^^. Apla- 
narity is defined as the smallest eigenvalue of the nor- 
malized laboratory-frame momentum tensor of the jets 
and the W boson. H!p2 = "^jp measures the event cen- 
trality, where Ht2 is the scalar sum of the transverse 
momenta of the jets excluding the leading jet and H\\ is 
the sum of the magnitudes of the momentum components 
parallel to the beamline of the jets, isolated lepton and 
neutrino. In this case the neutrino momentum parallel 
to the beam is estimated requiring that the mass com- 
puted from the measured lepton momentum, and un- 
known neutrino momentum parallel to the beam is equal 
to the W boson mass. If more than one solution is found, 
the one smallest in absolute value is used. The variable 

-^Tmin = — IS a mcasurc of the jet separation 

normalized by the transverse energy of the reconstructed 
W boson. Ai?™'" is the smallest distance in 77 — </) space 

between any two of the four leading jets, e}^^^'^'^ is the 
smaller of the two jet Ets. The transverse energy of the 

W boson is defined as E^ = p!^P*°" + |- These 
four variables are combined in a single discriminant vari- 
able -Dlb using the likelihood ratio procedure described 
in Ref. Q. 

For the analysis presented here, the low-bias discrim- 
inant -Dlb (= xi) was combined with a new variable 
called "pT-fraction" and the number of b tags to build 
a combined discriminant D. The pr-fraction, defined 

as X2 = (Etracks in jets PT)/(Eall tracks Pt), is the PT" 

weighted fraction of all tracks in the event that point 
to an energy deposit defining a jet (with jet pr >20 GeV 
with < 2.5). Only those tracks were considered that 
have a distance of closest approach of less than 1 cm 
along the beam direction with respect to at least one 
of the primary vertices in the event. In order to be in- 
cluded in the pt sum over tracks in a jet, a track was 
required to be within AR < 0.5 from the jet axis. This 
variable distinguishes clean events with nicely coUimated 
jets from events with broader jets and significant under- 
lying hadronic activity. Finally, 2:3 is the number of b 
tags. For each variable Xi, we use Monte Carlo simu- 
lation to determine the probability density functions Si, 
for tt signal, and bi, for iy-|-jets background. To a good 
approximation, these three variables Xi are uncorrelated, 
and the combined likelihood discriminant is derived as 

IliSi{xi)/bi{xi) 



D = 



UiSi{xi)/bi{xi 



1 



(3) 



thus combining event topology with a tracking-based jet 
shape and b tag information. This combined likelihood 
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discriminant offers a much better discrimination between 
tt and background than does the low-bias topological 
variable Dlb by itself, while maintaining its low level of 
correlation with the fitted top quark mass (and therefore 
with the jet energy scale). 

Figure [1] shows the distribution of the combined 
discriminant D obtained in the electron+jets and 
muon+jets channels. The distribution observed in data 
is compared to a model consisting of simulated tt and 
VF+jets events and the QCD multijets sample obtained 
from data fSec. lIVp. A likelihood fit is performed to de- 
termine the estimated fraction of events. The fit results 
are shown in Table HH In the fit, the ratio between the 
number of QCD and M^-|-jets events was kept fixed at a 
value based on the estimate used in Refs. ITL l23l. 



TABLE II: Composition of the 425 pb ^ data sample as 
determined by the likelihood fit. 





electron-fjets 


muon-fjets 


tt 


61.5 ± 7.9 


45.6 ± 7.5 


W+]ets 


35.6 ± 5.2 


63.0 ± 6.9 


QCD multijet 


18.9 ± 2.7 


5.4 ± 0.6 


Total observed 


116 


114 



from the constrained kinematic fit, which provides the 
sensitivity to the top quark mass and jet energy scale. 
To good approximation D is uncorrelated with xgt and 
with the jet energy scale. Thus the probabilities Pggn 
and Pbkg can be written as the product of a probability 
to observe a value D and a probability to observe xm , as 

Psgn (x; mt, JES) = Psgn (D) ■ P^gn (xRu mt, JES) (5) 

and 

Pbkg {x; JES) = Pbkg {D) ■ Pbkg (xfit; JES) (6) 

where D is calculated for a JES parameter equal to 1. 
The normalized probability distributions of the discrim- 
inant D for signal Pggn [D) and background Pbkg (P*) 
are assumed to be independent of JES and are ob- 
tained from Monte Carlo simulation as discussed in 
Sec. IVIIl They correspond to parameterized versions of 
the Monte Carlo templates shown in Fig. [1] The recon- 
struction of the signal and background probabilities for 
the mass information xgt is explained in Sec. lVIIl"Xl The 
mass information in the event Xfit consists of all fitted 
masses mi{JES), estimated uncertainties ai{JES), and 
goodnesses-of-fit x^{JES) obtained from the kinematic 
fit. 



VIII. THE IDEOGRAM METHOD 

To maximize the statistical information on the top 
quark mass extracted from the event sample, a likeli- 
hood to observe the event is calculated for each event as 
a function of the assumed top quark mass mj, the jet en- 
ergy scale parameter JES , and the fraction of ti events 
in the event sample, /top- The likelihood is composed of 
two terms, describing the hypotheses that the event is tt 
signal or background: 

£ovt (a;; mt, JES, /top) = /top • Psgn (a;; mt, JES) 

+ (1 - /top) • Pbkg (x; JES) . (4) 

Here, x denotes the full set of observables that charac- 
terizes the event, /top is the signal fraction of the event 
sample, and Pggn and Pbkg are the probabilities for tt 
and VF-|-jets production, respectively. The contribution 
from QCD multijet events is comparatively small and ex- 
pected to have a fitted mass shape very similar to that 
of VF-fjets events. Therefore no explicit QCD multijet 
term is included in the likelihood. The event observables 
X can be divided into two groups. One set is chosen to 
provide good separation between signal and background 
events while minimizing the correlation with the mass 
information in the event. These variables (topological 
variables and b tagging) are used to construct a low-bias 
combined discriminant Z?, as described in Sec. IVIIl The 
other event information used is the mass information Xj\t 



A. Calculation of signal and background 
probability 

The signal and background probabilities are calculated 
as a sum over all twenty-four possible jet/neutrino solu- 
tions. Without b tagging, the relative probability for each 
of the solutions i to be correct depends only on the xf for 
the corresponding fit and is proportional to exp(— ^x^). 
To further improve the separation between correct and 
incorrect jet assignments, b tagging is used. If one or 
more jets in the event are b tagged, an additional rela- 
tive weight Wbtag.i is assigned, representing the probabil- 
ity that the observed b tags are compatible with the jet 
assignment assumed for that particular jet permutation: 

wbtag,* = n c^) 

j = l,njct 

where pj can either be (1-e/), or (1-eh), depending 
on the assumed flavor of the jet (light or b) and whether 
or not that particular jet is tagged. For this purpose 
the jets from the hadronic W boson decay are always 
assumed to be light quark (u, d, s) jets. In the calibration 
of the analysis (see Sec. IVIlfE)) . however, the fraction of 
W cs decays and the higher tagging rate for c quark 
jets are taken into account. The tagging rates for light 
and b quark jets ei and £{, are used as parameterized 
functions of jet px and rj. The jet pt is based on the 
reconstructed jet energy for JES = 1, consistent with 
the jet energy scale for which the tagging rate functions 
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FIG. 1: Combined likelihood discriminant D in data and MC simulation in the electron+jets channel (left) and muon+jets 
channel (right). The tt, W^+jets and multijet contributions are normalized according to the fitted fractions. 



are derived from data [17]. Thus, the weight assigned to 
each jet combination becomes 



= CXp(--X,) • Wbtag. 



(8) 



The mass-dependent signal probability in Eq. [5] is calcu- 
lated as 



-Psgn {xau nit, JES) 



E 



Wi 



G{mi, to', ffi) ■ BW(m', mt)dm! 



(1 /correct) ' ^wronni^ii ''^t) 



(9) 



The signal term consists of two parts: one part de- 
scribes the compatibility of the solution with a certain 
value of the top quark mass, assuming that it is the 
correct solution. It takes into account the estimated 
mass resolution (T; for each jet permutation. The sec- 
ond part of the signal term describes the expected shape 
of the mass spectrum for the "wrong" jet assignments, 
which also depends on the top quark mass. The "cor- 
rect" solution part is given by a convolution of a Gaus- 
sian resolution function G{mi,m' ,ai) and a relativis- 
tic Breit-Wigner BW(TO',TOt). The Gaussian function 
describes the experimental resolution. The relativistic 
Breit-Wigner represents the expected distribution of the 
average invariant mass m' of the top and anti-top quark 
in the event for a given top quark mass toj. The width 
of the Breit-Wigner is set to the standard model value 
of the top decay width Q. The "wrong" permutation 
signal shape S^^^ng{mi,'mt) is obtained from MC simu- 
lation using a procedure described in Sec. IVIII Bl These 
two terms are assigned relative weights depending on 
/corTcct' which represents the relative probability that the 
weight is assigned to the correct jet permutation. For 



well-reconstructed events with exactly 4 jets, this proba- 
bility is approximately 39% if b tagging is not used. For 
4-jet events with 0, 1, or > 2 tagged jets, the values 
/correct— 0-45, /4,rrcct=0.55, and /^Q^^g(.t=0.65 are used. 
For 5-jet events smaller fractions are used: 0.15, 0.30, 
and 0.40 for events with 0, 1, or > 2 6 tagged jets re- 
spectively. Ensemble tests (see Sec. IVIII Ep confirm that 
these values result in a pull width for the mass close to 
unity for the different tagging multiplicities. 
The background term in Eq. [S]is calculated as: 



-Pbkg (x^t; JES) = ^Wi- BG{mi) 



(10) 



1=1 



where the background shape BG{m) is the shape of 
the fitted mass spectrum for V7-|-jets events. To obtain 
BG{m), the kinematic fit (with JES equal to unity) is ap- 
plied to simulated W-|-jets events and the fitted masses 
mi for all possible jet/neutrino solutions i are plotted. 
All entries are weighted according to the permutation 
weight Wi defined in Eq. [5] The shapes of BG{m) used 
in the analysis are shown in Fig. O 

The Breit-Wigner and "wrong" permutation signal 
shape are normalized to unity within the integration in- 
terval of TOmin = 100 GeV to TOmax — 300 GcV. This 
interval is chosen to be large enough so as not to bias the 
mass in the region of interest. 

The normalization of the background shape _BG(to) is 
chosen such that the fitted signal fraction /top reproduces 
the true tt fraction in ensemble tests (see Sec. IVIII E[) 
containing ti and M^-|-jets events. The mass fit tends to 
underestimate /top, due to the presence of ti events that 
are misreconstructed or affected by energetic gluon radi- 
ation and resemble W^-|-jets events in the fact that their 
topology does not conform to the ti hypothesis in the 
kinematic fit. A constant normalization factor of 1.15 is 
found to reduce the offset in /top to less than 1% both 
in the electron-|-jets and the muon-|-jets channel. The 
jet energy scale parameter is varied before performing 
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FIG. 2: Histograms showing the background shape from a weighted sum (see text) of all twenty-four masses from each event 
from the W+jets background sample (points with error bars), for the electron+jets channel (left) and muon+jets (right). The 
histograms show the shapes that are used in the likelihood. To reduce statistical fluctuations, the shapes are calculated as the 
average value in a sliding window of ±5 GeV around each fitted mass. 



the constrained fit by scaling all jet energies by a con- 
stant factor. The event likehhoods are recalculated for 
each different value of the JES parameter. Since the con- 
strained fit uses a W boson mass constraint, the 
the fit will be best when the invariant mass of the jets 
from the hadronically decaying W boson is closest (on 
average) to the known W boson mass. Additional sensi- 
tivity to the jet energy scale comes from the shape of the 
fitted mass distribution in background events. For the 
proper jet energy scale the spectrum will agree best with 
the background shape included in the background term 
in the likelihood. 



B. Determination of the wrong-permutation signal 
shape 

The convolution of Gaussian detector resolution and a 
Breit-Wigner, used in the signal term of the likelihood, 
implicitly assumes that the correct jet assignment is cho- 
sen. To describe the contribution from wrong jet assign- 
ments, a separate term is added to the signal part of the 
likelihood. To obtain the fitted mass spectrum of the 
wrong permutation signal, samples of parton-matched 
tt events are used in which all quarks are matched to 
jets. The fitted mass spectrum is plotted including all 
jet permutations except the correct solution (excluding 
both neutrino solutions corresponding to the correct jet 
permutation). Each entry is weighted according to the 
permutation weight assigned in the Ideogram likelihood. 
Samples of different generated top quark masses are used. 
For each mass, the weighted sum of wrong solutions is fit- 
ted with a double Gaussian. The fitted parameters for 
correct solutions and for the wrong permutation signal 
show a linear behavior as a function of the top quark 
mass. The fitted parameters are given in Table Hill Since 
the permutation weights change when b tagging is in- 
cluded, this exercise is repeated for events with tags, 1 



tag, and 2 or more tags. 

The linear fits are used to construct a 2-dimensional 
wrong-permutation signal shape as a function of the fit 
mass and generated top quark mass S^^^„g{mi, rrit). For 
each value of the generated top quark mass, the shape as 
a function of fitted mass is described as the sum of two 
Gaussians. The resulting parameterizations are displayed 
as the wrong-permutation shapes in Fig. [3] and ID Also 
shown are the shapes of the correct jet assignments, de- 
termined in a similar fashion from parton-matched events 
using a single Gaussian. A linear dependence of the pa- 
rameters is found as a function of generated top quark 
mass. The sum of the correct solutions and wrong so- 
lutions is compared to a weighted histogram of all fitted 
masses in tt simulation. The parametrized functions give 
an adequate description of the overall (wrong -I- correct) 
signal shape. In Fig. [31 the corresponding distributions 
are shown for events with 0, 1, or 2 tags. It is clearly 
visible how the fraction of the weight given to the cor- 
rect solution improves when including b tag information 
in the permutation weights. In Fig. 31 the nine distribu- 
tions are shown for generated top quark masses ranging 
from 150 GeV to 200 GeV. 



C. Determination of JES offset correction 

The likelihood fit relies on the invariant mass of the 
hadronically decaying W boson in the ti events to set 
the jet energy scale. It is designed to give an unbiased 
fit of the JES parameter in well-reconstructed ti events 
when the correct jet assignment is used. However, in a 
significant fraction of the events, the jets that are pre- 
sumed to originate from the W boson may not really 
come from a W boson. Such cases include events other 
than tt, as well as tt events that are mis-reconstructed. 
In the presence of such events we can expect an offset 
in the fitted JES parameter. The slope of the JES call- 
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TABLE III: Parameters used to describe the background shapes (arbitrary normalization). For each case, the shape is described 
by the sum of two Gaussians G'(mflt) = a ■ exp [— (p — mat)^ /2a'^] , where the three parameters a, /i, and a evolve linearly as a 
function of the generated top quark mass mt as pO + pl • {nit — 175 GeV). 

tags 1 tag >2 tags 

Gaussian 1 Gaussian 2 Gaussian 1 Gaussian 2 Gaussian 1 Gaussian 2 

Parameter pO pi pO pi pO pi pO pi pO pi pO pi 

a 284.9 -1.722 51.72 -0.4199 267.5 -1.0700 68.08 -0.7129 235.5 -0.1662 75.86 -0.0415 

fi 161.7 0.7383 223.1 1.242 162.6 0.7859 220.1 1.400 166.2 0.6416 229.4 0.7454 

cr 23.55 0.2392 22.94 -0.2528 23.27 0.2737 23.97 -0.4551 25.80 0.1165 21.78 -0.2828 
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FIG. 3: Prediction of the shapes of the fitted mass distribu- 
tion for the wrong and the correct permutations (hatched) 
and the sum of the two (black line) using the fitted param- 
eters shown in Table IIIII The sum of the two is compared 
to the simulated data containing a weighted sum of all so- 
lutions (correct and wrong), for the default jet energy scale 
and a generated top quark mass of 175 GeV for events with 
(upper), 1 (middle), or more than 1 (lower) b tags. 



bration curve (fitted JES parameter as a function of the 
"true" JES) may also differ from unity. 

Using the MC calibration procedure described in 
Sec. IVIII E[ we find that the presence of wrong jet as- 
signments and background events causes an offset of sev- 
eral percent in the fitted JES parameter. A breakdown 
of the different contributions to the JES offset and slope 
is shown in Table IIVI 



TABLE IV: The JES calibration slope and offset for different 
event samples are shown. The offset increases and the cali- 
bration slope becomes smaller when mis-reconstructed signal 
events or background events are added. The offset correc- 
tion at the likelihood level (see text) fixes the JES offset but 
further reduces the JES calibration slope. 





JES slope JES offset 


Smt (*) 


parton-matched tt only 


0.96 


+0.026 




ti only 


0.88 


+0.050 




all events 


0.80 


+0.076 


4.30 GeV 


all, 50% offset correction 
all, 100% offset correction 


0.72 
0.63 


+0.036 
+0.000 


4.10 GeV 
4.01 GeV 



(*) expected mass uncertainty after full calibration 



The JES offset and slope turn out to be independent 
of the generated top quark mass (see Fig. [5]). Therefore 
we apply a straightforward mass-independent correction. 
A normalization factor fjEsiJES, /top) = exp(a • JES) is 
introduced which corrects for the offset without changing 
the statistical uncertainty estimated from the likelihood 
(in case the final sample likelihood is Gaussian): 

£™7(m^, Ji?5, /top) = 

fjEsiJES, /top) • £ovt(mf , JES, /top). (11) 

Since background events on average cause a larger bias 
than signal events, a is defined to be dependent on the 
measured signal fraction /top: a = 2.63 + 0.56(1 — /top)- 
The value of the correction constant is tuned using MC 
simulation to give an unbiased measurement of the JES 
at the reference scale JES — 1. As shown in Table IIVI 
the application of this offset correction removes the JES 
offset, but it further reduces the JES calibration slope. 
Table IIVI also shows that after full calibration (described 
in the next section), the expected statistical uncertainty 
on the top quark mass improves slightly when applying 
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FIG. 4: Same as Fig. [S] for difFerent values of the generated top quark mass, combining all events irrespective of the number 
of b tags. 



the corrections. For illustrative purposes we also include 
a 50% offset correction in the table, where 0.5 • a is used 
instead of a. 

The correction described above ensures that the fit is 
well-behaved and that, for values of the JES parameter 
near 1, the fit results will stay well within the range for 
which the {JES, nit) likelihood is calculated. It does not, 
however, provide a full calibration of the analysis, which 
is described in Sec. IVIII El 



D. Combined likelihood fit 

Since each event is independent, the combined likeli- 
hood for the entire sample is calculated as the product 
of the single event likelihood curves: 

£samp(mt, JES, /top) - n ^c°tj{mt, JES, /top). (12) 

j 

This likelihood is maximized with respect to the top 
quark mass rrit, the jet energy scale parameter JES, and 
the estimated fraction of signal in the sample /top. 



E. Calibration using Monte Carlo simulation 

The analysis is calibrated using Monte Carlo simula- 
tion. Both the bias on the measured mass and the cor- 
rectness of the estimated statistical uncertainty are stud- 
ied using ensemble tests, in which many simulated ex- 
periments (pseudoexperiments) are created, each match- 
ing the size of the observed data sample. Thousands 
of pseudoexperiments are constructed, combining tt and 
W^-|-jets events from MC simulation. The event fractions 
for tt and VF-|-jets are allowed to fluctuate according to 
binomial statistics around the estimated fractions in the 
actual data sample. The fractions used are those listed 
in Table |lll In the pseudoexperiments, the QCD mul- 
tijets contribution is replaced by VF-|-jets events. This 
deviation in QCD multijet fraction is treated as a sys- 
tematic uncertainty (see Sec. [X]). The total sample size 
is fixed to the observed number of events in data (116 in 
electron-|-jets and 114 in muon-hjets). To make optimal 
use of the available MC statistics, standard resampling 
techniques are used, allowing for the multiple use of MC 
events when constructing the pseudoexperiments [2^ . 
For every pseudoexperiment the mass is fitted and the 
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FIG. 5: The mean fitted JES and pull width as a function of the generated top quark mass Mgen for a "true" JES of 0.97 
(left), 1.00 (middle), and 1.03 (right), for the lepton-l-jets channel (e -|- At combined). The fitted JES is stable as a function of 
generated top quark mass. 
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FIG. 6: The difference between the mean fitted mass Mfit and the generated top quark mass Mgen as a function of the generated 
top quark mass for a "true" JES of 0.97 (left), 1.00 (middle), and 1.03 (right), for the lepton-|-jets channel (e -|- /i combined). 
At a generated mass of 175 GeV, the mass bias changes by 1 GeV when the true JES is varied by ±3%. 
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deviation of this mass from the mean of all pseudoexper- 
iment masses is divided by the fitted uncertainty. This 
quantity is referred to as the "pull." The pull distribu- 
tion for all pseudoexperiments is fitted with a Gaussian 
to extract the width, which we call the "pull width." The 
corresponding pull and pull width for the fitted JES are 
also determined. 

Figure [5] shows how the mean fitted JES and its pull 
width behave as a function of the top quark mass for dif- 
ferent values of the true jet energy scale. The fitted JES 
parameter is independent of the top quark mass over the 
full range considered. The plots also show that the fitted 
JES changes linearly as a function of the true JES with 
a slope of 0.63 (see discussion in Sec. I VIII C|) . Figure [6] 
shows the change in the fitted top quark mass and the 
width of the pull as a function of the generated top quark 
mass for different values of the true JES . Using these 
plots a full two-dimensional calibration is performed, de- 
scribing the fitted JES and top quark mass as a function 
of the "true" JES and top quark mass generated in the 
MC simulation. The estimated statistical uncertainties 
are corrected for the width of the pull and error propa- 
gation is used to take into account the effect of the two- 
dimensional calibration including the correlation between 
the JES parameter and the offset in measured mass. 



F. Alternative JES fitting strategies 

Including the uniform JES parameter as a free param- 
eter in the fit reduces the systematic uncertainty due to 
the jet energy scale, at the cost of a larger statistical 
uncertainty. As a comparison, in Fig. [7] the expected sta- 
tistical uncertainties on the top quark mass are shown 
for three different fitting scenarios. When fixing the JES 
parameter in the fit to 1, the statistical uncertainty is 
smallest: 2.59 GeV at a generated top quark mass of 
175 GeV. When allowing the JES parameter to float 
freely in the fit, without correcting for the JES slope 
(=0.63) in the calibration, part of the sensitivity to the 
overall JES scale is reduced and absorbed in the statis- 
tical uncertainty, leading to an expected statistical un- 
certainty of 3.34 GeV at a top quark mass of 175 GeV. 
Fully calibrating the analysis as a function of fitted mass 
and JES (the default approach), allows an unbiased top 
quark mass measurement for any value of the "true" JES , 
at the cost of a larger statistical uncertainty: 4.01 GeV 
at a generated top quark mass of 175 GeV. 

In order to be consistent with the approach used by 
the Matrix Element analysis 0], thus facilitating a com- 
bination of results, and to minimize the dependence on 
the external JES constraint from jet+photon studies, the 
third scenario is presented here as the main analysis re- 
sult, applying the full calibration as a function of fit- 
ted top quark mass and JES. Results using the other 
two JES fitting strategies are quoted as a cross-check in 
SecKH 



IX. RESULTS WITH DATA 

The overall likelihood curves obtained for data are 
shown in Fig. [S] The 2D likelihoods show the actual 
likelihood values in bins of 1 GeV in mass and 3% in 
JES. The jagged appearance of the ellipses is caused by 
the large bin size in the JES direction. To extract the 
mass and statistical error, a Gaussian fit is applied to the 
three bins closest to the minimum in the one-dimensional 
negative log likelihood curves. The fitted values are cor- 
rected according to the calibration derived in Sec. lVIIIEl 
The measured top quark mass is: 

mt = 173.7 ± 4.4 (stat -I- JES) GeV 
with 

JES = 0.989 ±0.029 (stat). 

All uncertainties shown are statistical. The fitted ti sig- 
nal fraction is /top = 0.453 ±0.032. If the JES parameter 
is kept fixed to 1 in the fit, the estimated statistical un- 
certainty is 2.93 GeV. Hence the 4.43 GeV (stat±JES) 
uncertainty of the 2D fit can be interpreted as a combi- 
nation of an intrinsic mass uncertainty of 2.93 GeV (stat) 
and an additional uncertainty of 3.32 GeV (JES) due to 
fitting the JES parameter. As shown in Fig. [51 the ob- 
served statistical uncertainties are slightly larger than the 
average uncertainties expected from Monte Carlo ensem- 
ble tests, but they fall well within the distribution. The 
fitted JES of 0.989±0.029 is in good agreement with the 
reference scale 1 (or 5=0), corresponding to the hypoth- 
esis that after all jet corrections the JES in data and MC 
are the same. 

One can also compare the in-situ fitted JES parame- 
ter with the scale obtained in jet+photon studies. When 
correcting all jets in MC events for the jet-p^ depen- 
dent difference between data and MC, 5, and redoing 
the ensemble tests in MC simulation, the mean fitted 
JES is 0.962^Q Q23, where the uncertainties correspond 
to the combined statistical and systematic bounds from 
the jet+photon studies. This is consistent with the value 
of 0.989 ± 0.029 measured in situ. 



X. SYSTEMATIC UNCERTAINTIES 

The calibration of the analysis relies on Monte Carlo 
simulation. Therefore any discrepancy between the 
Monte Carlo simulation and the data may lead to a bias 
and thus to a systematic shift in the measured top quark 
mass. In this section we describe the aspects of the sim- 
ulation which may not accurately represent the data and 
evaluate the possible effect on the mass measurement. 
To determine the impact of each uncertainty, we perform 
ensemble tests using a pool of simulated events that are 
modified according to the uncertainty in question. The 
shift in the mean fitted ensemble mass compared to the 
default value gives the size of the corresponding system- 
atic uncertainty. 
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FIG. 7: The expected statistical uncertainty from ensemble tests is shown as a function of the generated top quark mass 
for three scenarios: with the JES parameter fixed to 1 (left), allowing the JES parameter to float freely in the fit, but only 
calibrating the mass fit for a true JES — 1 (middle), allowing the JES parameter to float freely in the fit and applying the 
full calibration as a function of true top quark mass and true JES (right). In each plot the width of the band indicates the 
estimated uncertainty. 
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FIG. 8: Overall likelihood curves for the events observed in data, in the electron+jets channel (left), muon+jets (middle), 
and both channels combined (right). The top plots show the full 2-dimensional likelihood as a function of the jet energy scale 
parameter (JES) and top quark mass. Each contour, n, corresponds to a difference in likelihood of Aln(/I) = —n?/2 with 
respect to the maximum likelihood. The fitted value of the JES parameter as a function of the top quark mass is plotted as the 
gray line superimposed on the 2D likelihoods. The bottom plots show the likelihood as a function of the top quark mass along 
the gray line from the upper plots. The fitted values from these distributions have to be corrected for the calibration from MC 
simulation to obtain the final results. 
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FIG. 9: Distribution of the estimated statistical uncertainty on the top quark mass measurement (left) and JES measurement 
(right) for the fully calibrated analysis, in the combined lepton+jets channel. The values observed in data are indicated by the 



The total systematic uncertainty on the top quark 
mass measurement is obtained by adding all contribu- 
tions in quadrature. The following sources of systematics 
are considered (also see Table W} : 

• Jet energy scale px dependence: The inclusion 
of a uniform jet energy scale parameter JES as a 
free parameter in the mass fitting and calibration 
procedure ensures that a relative difference in over- 
all jet energy scale between data and Monte Carlo 
is corrected for. The corresponding uncertainty is 
included in the quoted statistical (stat -|- JES) un- 
certainty. Any residual discrepancy between data 
and Monte Carlo jet energy calibration that can- 
not be described by a uniform scale factor may lead 
to an additional systematic error on the top quark 
mass. The largest additional effect is expected from 
the uncertainty in the jet-px dependence The 
size of the impact of a possible jet-px dependent 
shape is estimated by scaling the energies of all 

jets in the MC with a factor (1 + 0.02^^^^^^^), 

where p^^* is the default reconstructed jet pt- The 
value of 0.02 is suggested by the jet+photon stud- 
ies. The mass obtained with the modified pseudo- 
experiments is compared to the default result and 
the shift of 0.45 GeV is quoted as a systematic un- 
certainty. 

• Jet reconstruction efficiency and resolution: 

In addition to uncertainties on the reconstructed jet 
energies, differences between data and the Monte 
Carlo simulation in the jet reconstruction efficiency 
and jet energy resolution may lead to a mass bias. 
Both efficiency and resolution are varied as a func- 
tion of jet Px and rapidity within estimated un- 
certainties. No significant effect is observed, with 
an estimated statistical precision of 0.15 GeV. For 
both effects combined, a systematic uncertainty of 
0.22 GeV is quoted. 

• 6-fragmentation: While the overall jet energy 
scale uncertainty is included in the statistical un- 



certainty from the fit, differences between data and 
Monte Carlo in the ratio of 6-jet and light-jet en- 
ergy scale could still affect the measurement. One 
possible source for such differences could be the 
description of 6-jet fragmentation in the simula- 
tion. To estimate the uncertainty from this source 
we used samples of simulated tt events with dif- 
ferent fragmentation models for b jets. The de- 
fault Bowler [2^ scheme with rti=1.0 is replaced 
with rb=0.69 or with Peterson [26| fragmentation 
with £b=0. 00191. These parameter values were 
obtained by tuning pythia simulation to LEP 
data [13, H^, [2^]. The size of the variation in rt 
corresponds to a larger shift in mean scaled energy 
(xb) of b hadrons than the uncertainties reported 
in ^3, m, 29, 30]. The comparison between the 
Bowler and Peterson scheme addresses the uncer- 
tainty on the shape of the xb distribution. Ensem- 
ble tests are repeated using events from each of the 
three simulations. The absolute values of the de- 
viations in top quark mass results with respect to 
the standard sample are added in quadrature and 
quoted as a symmetric uncertainty of 1.3 GeV. 

• b-jet energy response: Uncertainties in the sim- 
ulation of the ratio between the calorimeter re- 
sponse to hadronic showers and electromagnetic 
showers (h/e ratio) may lead to additional differ- 
ences in the &/light jet energy scale ratio between 
data and simulation. The possible size of the ef- 
fect is studied in simulation, combining the uncer- 
tainty in the modeling of h/e calorimeter response 
ratio with the difference in particle content between 
light quark and b jets. An estimated uncertainty of 
1.4% on the 6/light jet energy scale ratio is found. 
Ensemble tests show that this corresponds to a 
systematic uncertainty on the top quark mass of 
1.15 GeV. 

• b-tagging: The b tagging rates for b jets, c jets, 
and light-quark jets are varied within the uncer- 
tainties known from the data, and the resulting 
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TABLE V: Summary of systematic uncertainties. 

Source of uncertainty Size of the effect (GeV) 

Jet energy scale {pr dependence) 0.45 

Jet ID efficiency and resolution 0.22 

6 fragmentation 1.30 

6 response (h/e) 1.15 

b tagging 0.29 

Trigger uncertainty +0.61 — 0.28 

Signal modeling 0.73 

Signal fraction (stat+sys) 0.12 

Background modeling 0.20 

Multijet background 0.28 

MC calibration 0.25 

PDF uncertainty 0.023 

Total systematic uncertainty +2.10 — 2.04 



variations are propagated to the final mass results. 
Uncertainties in the heavy flavor composition of the 
background are also considered. The combined ef- 
fect is 0.29 GeV. 

• Trigger: The trigger efficiencies in the Monte 
Carlo simulation are varied by their uncertainties 
estimated from data. The resulting variations in 
fitted mass are summed in quadrature, leading to a 
combined trigger uncertainty of +0.61 — 0.28 GeV. 

• Signal modeling: The main uncertainty in the 
modeling of tt events is related to the radiation of 
gluons in the production or decay of the tt system. 
A difference in the description of hard gluon radi- 
ation could affect the transverse momentum spec- 
trum of the tt system or, for example, change the 
rate of confusion between jets from the hadronically 
decaying W boson and initial state gluons, which 
could affect the reconstructed top quark mass. To 
assess the uncertainty related to the modeling of 
high energy gluons, the difference is studied be- 
tween the default signal simulation and a dedicated 
tt+jet simulation in which an energetic parton is 
produced in addition to the system in the produc- 
tion process simulated by ALPGEN. ft is estimated 
that in the class of events that pass the full event 
selection, the fraction of simulated events with such 
an energetic gluon disagrees with the data by less 
than 35%. Pseudoexperiments are made with the 
usual sample composition, but replacing the default 
tt events with the events from the dedicated tt+jet 
simulation. 35% of the observed shift in the fitted 
mass corresponds to 0.73 GeV, which is assigned as 
a systematic uncertainty. 

• Signal fraction: Since the tt fraction /top is fit- 
ted together with the top quark mass and the JES 
parameter, the mass measurement is affected by 
the uncertainty on the signal fraction in the data 



sample. We estimate two sources of systematic un- 
certainty: a variation of the signal fraction in the 
ensemble test used to calibrate the method and the 
effect of a possible systematic offset in the fitted sig- 
nal fraction with respect to the true signal fraction 
internally in the mass fit. 

We take the 7% relative statistical uncertainty of 
the signal fraction found by combining the /i+jets 
and e+jets numbers shown in Table |TT1 We add 
in quadrature an estimated relative systematic un- 
certainty of 11% estimated from the cross section 
measurements [l3l ■ New ensemble tests for the cal- 
ibration procedure, are performed with the mean 
of the Poisson distribution (for the signal fraction) 
shifted by (11 © 7)%. Following this procedure the 
combined fit will still correctly fit the different sig- 
nal fraction and compensate for the effect. This 
does not take into account the effect of a possi- 
ble systematic discrepancy between the data and 
the Monte Carlo model of signal and background, 
which could lead to a systematic bias in the fitted 
signal fraction. To evaluate this additional system- 
atic, the mass fit is forced to systematically over- 
estimate or underestimate the tt fraction by 11% 
(with respect to the value preferred by the likeli- 
hood fit) , and the shift in fitted mass is quoted as a 
systematic uncertainty. The combined uncertainty, 
adding the above two contributions in quadrature, 
is 0.12 GeV. 

• Background modeling: The sensitivity of the 
measurement to the choice of background model 
is studied by comparing two enlarged pseudoex- 
periments in which the background simulation is 
changed. One sample is based on the standard 
W^+jets simulation using a factorization scale of 
Q"^ = + X^jIPt)^ while in the other pseudoex- 
periment a sample of W^+jets events is used that 
are generated with a different factorization scale 
of Q'^ = (Px)"^- The observed difference in fitted 
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mass is 0.20 GeV which is assigned as a systematic 
uncertainty. 

• QCD multijet background: In the cahbration 
procedure, the M^+jets simulation is used to model 
the small multijet background in the selected data 
sample. To study the systematic uncertainty due 
to this approximation, we selected a dedicated 
multijet-enriched sample of events from data by in- 
verting the lepton isolation cut in the event selec- 
tion. The calibration of the method is carried out 
with pseudoexperiments in which these events are 
used to model the multijet background, according 
to the fractions given in Table [TTl The observed 
shift is 0.28 GeV, which is quoted as a systematic 
uncertainty. 

• MC calibration: The statistical uncertainty on 
the calibration curves shown in Fig.[n]is propagated 
through the analysis and yields a systematic uncer- 
tainty on the result of 0.25 GeV. 

• Uncertainty due to the parton distribution 
functions (PDF): The Ideogram analysis mea- 
sures the top quark mass directly from the invari- 
ant mass of the tt decay products without making 
specific assumptions regarding the production pro- 
cess. Nevertheless, the calibration of the analysis 
relies on Monte Carlo simulation in which a certain 
PDF set was used (CTEQ5L [31]). It is conceiv- 
able that a different choice of PDFs would lead to a 
slightly different calibration. To study the system- 
atic uncertainty on the top quark mass due to the 
precise PDF description, several PDF uncertainties 
are considered. PDF variations provided with the 
next-to-leading-order PDF set CTEQ6M [s^ are 
compared to the default CTEQ6M. The difference 
between CTEQ5L and MRST leading order PDFs 
is taken as a separate contribution. Also the ef- 
fect of a variation in as is evaluated. In all cases 
a large pseudoexperiment composed of events gen- 
erated with CTEQ5L is reweighted so that distri- 
butions corresponding to the desired PDF set are 
obtained. The difference between weighted and 
unweighted pseudoexperiments is then quoted as 
systematic uncertainty, and all individual uncer- 
tainties are added in quadrature. The resulting 
combined uncertainty is found to be very small: 
±0.02 GeV. 



XI. CROSS-CHECK USING AN EXTERNAL 
JES CONSTRAINT 

As a cross-check, the analysis is repeated using the two 
alternative JES fitting strategies discussed in Sec. lVIIIFl 
Fixing the JES parameter in the fit and relying fully on 
the external JES constraint from jet+photon studies, the 



top quark mass is measured to be: 

mt = 175.8 ±2.9 (stat)+2:7 (JES) GeV, 

quoting only the statistical uncertainty (stat) and the 
systematic uncertainty due to the jet energy scale (JES). 
In the other alternative approach, the JES parameter 
is allowed to float freely in the fit but no calibration of 
the JES slope is applied. Again, the external JES con- 
straint from jet+photon studies is required to set the jet 
energy scale and the remaining JES systematics. Effec- 
tively this approach combines in-situ with external JES 
information, leading to the following result: 

mt = 173.9 ±3.6 (stat)^^ ^ (JES) GeV. 

Comparing the last (most precise) cross-check with the 
main result, one can conclude that omitting the external 
JES constraint and relying fully on the in-situ informa- 
tion changes the central result only by 0.2 GeV. The 2 
GeV difference between the first cross-check and the main 
result correlates very well with the 1.1% difference in JES 
value between the default Monte Carlo scale and in-situ 
JES measurements. This difference is fully covered by 
the quoted uncertainties. 

XII. CONCLUSION 

The Ideogram method has been used for the first time 
to measure the top quark mass in tt events with the 
£+jets topology. This technique employs a kinematic 
fit to extract mass information from the events, while 
improving the statistical sensitivity by constructing an 
analytic likelihood for every event taking into account 
all jet permutations and the possibility that the event is 
background. Lifetime-based identification of b jets is em- 
ployed to enhance the separation between ti signal and 
background and to improve the assignment of the ob- 
served jets to the partons in the ti hypothesis. To reduce 
the systematic uncertainty due to the jet energy scale 
calibration, an overall scale factor JES for the energy of 
the reconstructed jets is a free parameter in the fit deter- 
mined simultaneously with the top quark mass and the 
signal fraction. 

From a DO Run II data sample of approximately 
425 pb^^, 116 events are selected in the electron+jets 
channel and 114 in the muon+jets channel. The top 
quark mass is measured to be 

mt = 173.7 ± 4.4 (stat + JES);2:o (^yst) GeV 

with a fitted JES scaling factor: 

JES = 0.989 ± 0.029 (stat only), 

which is consistent with the reference jet energy scale 
(=1.0) and with the results from the jet+photon calibra- 
tion (« 0.9621q q23). The mass result is in good agree- 
ment with the Matrix Element measurement using the 
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same data set [Tj and with other recent top quark mass 
measurements [J, . 
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